UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRARY ROUTINE V 6 - 202, 

M l l q 



TITLE 

TYPE 

NUMBER OF WORDS 

TEMPORARY STORAGE 

PRESET PARAMETER 

DURATION 
ACCURACY 



■where 



Associated Legendre Functions 

Closed. 

80 (71 plus the 9 of Program R 1 -which is 

attached . ) 

0-10 inclusive and (k + 2) ( £ + l) words at 

n n 

S3 

S 3 contains the address at -which the table of 

functions IT: (x), K- (x), KT (x)* ... K n (x) is 

to be placed 

[25 + 24k + 12 k ]ms , Approximately (this 

expression holds if k '> 2) 

The maximum errors that can occur due to round 

off are: 



__ m _ 
e* <« error in j£ ) . m2" 59 + I ff) x J ^ 



j=0 



-m+j 



,0 



,-59 



k-2 



€^ (= error in l£ ) = 2"^ E (2X) 1 

i=0 

(Note that e°, e° < 2**^) 



DESCRIPTION AND ENTRY 



In general the errors are reasonably less than 
these maxima. See lis cuss ion of Error below « 

Enter with -*985 < x < -2*" 19 , 2~ 19 < x < .99 in 
the accumulator and the program parameter: 



< k < 12 



p 






50 
50 


k F 

" n 

pF 


p 


+ 


1 


26 


qF 



\z 



where k is the maximum order of the functions 
desired, q. is the adilress of this routine. KT (x) will be placed in the 



- 2 - 

location [l/2 (k + l) (k) + m] S3, i,e. K° in S3, K^ in 1S3, K^ in 2S3, 
£j in 3S3, k£ in 4S3, K^ in 5S3, • •, E^ n in [1/2 (k n + 2)(k n + l)] S3 
If k > 13 the functions KT will exceed scale. If k is 0, Ky K, , 
and K7 will be computed. 

NOTE This program contains the library routine R1- 116 

Square Root as an appended part. It is hence available for use in other 
connections and would not need to be placed a second time in the machine. It 
can be entered in the usual fashion at word 71L of this routine. 

DISCUSSION OF METHOD ™ „ , . T Jii / N , , ' ., . . . 
^iu^uuivu ^ i* f jt iie fyjjctions KT (x) computed by this routine are 

defined as 



(x) = _i_ ff ( X) = Sai. JtsilL. ^ (x) (d 

2/T~ * V Bjt (k+m)l * 
where P, (x) are the standardly defined Associated Legendre Polynomials* p™ (x) 
are the Normalized .Associated Legendre Polynomials so that K, satisfy: 

#*)£(x)a*= ^- V (2) 



r 



For all Q. < m < k and all < k < 12, the functions |K? (x) | < | and are 
hence in scale. 

The routine operates in three parts: 
(a) The unnormalized Legendre Polynomials P - P are first computed 
using: 

1/2 P Q (x) = 1/2 

1/2 P (x) = 1/2 x 

1/2 P k (x) = x[l/2 P k . 1 (x)] + ^- £[1/2 P^(x)]-l/2 P^W) 

(3) 

and inserting l/2 P (x) in the proper place in the table at S3» 
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(b) Each l/2 P k (x) is multiplied "by 2 J— ^ to form 



<m- 2 I^ ri/ 2 p k (x)] w 



(c) The recursion formula 



^—2^ ™ (x) = 1 (2k+l)(b*.l) ^-1 (x) . /^2±1 x j^-1 (x) (5) 

V (2k-l)(k+m) V k + m 

is then used to compute the K? (x) for M - and fill in the table as S3* 

DISCUSSION OF ERROR This program suf f er s from the very great defect 
common to all programs which generate a sequence of functions using recursion 
relations 9 namely the accumulation and magnification of error at each step. 
From equations (3) and (5) it can be seen that the error in K increases with k 
and that the error in K/ for fixed k increases with m. The latter increase 
is much more pronounced tha n the fo rmer* In addition the error will be much 
worse for | x | ~ 1 due toll - x in equation (5)* In fact for x < - .985 or 
x > .99, the K? lose so much significance for k and m large that division 



... -19 
hang-ups will occur, (if | x | < 2 a division hang -up will occur in the 



attempt to compute 1 - x ). These possibilities dictate that this routine be 

used with a great measure of care, skill, and caution. Representative 

K?J calculated by this routine for representative values of x have been checked 

against the 10 place tables of Associated Legendre Function of Zaki, Mtirsi, 

-10 
Fouad I University, 19^1 • These errors are given in Table I in units of 10 

Note the extremely poor results for x near 1 and for k, m large. One saving 

feature in this respect is that K? (for m \ 0) -» as |x | -» 1. Hence in the 

range near 1 in which this program is inapplicable, it may be permissible in 

some problems to set KT = for x in this troublesome range. For a fixed 

k this becomes a better approximation as m increases which arejust the cases 

(i.e. m increasing for fixed k) in which this routine becomes increasingly 

bad . 
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TABLE.- I: ERROR IN TBS FUNCTIONS E^ (in units of 10" 10 ) 


FUNCTION 


X = 


JC = -3 


X = .6- 


x = ,9 


X = .95 


X = *98 


x = .985 


4 


< 1 


< 1 


3 


9 


1 


7 


2 


4 




3 


1 


10 


9 


5 


8 


4 


<^ 


< l 


l 


*3 


6 


'; 2 


2 


4 


<1 


<< l 


< l 


6 


< 1 


1 


7 


4 




< 1 


< l 


< 1 


.1 


A 


2 


3 
K 3 


< i 


< 1 


< l 


•7 


2 


9 


11 


r 1 




I 


< 1 


< 1 


4 


4 


I 


^4 


< i 


.3 


< 1 


« 1 


2 


10 


1 


4 


< i 


< 1 


< l 


2 


9 


50 


62 


4 


< i 


4 


< l 


2 


6 


5 


2 


4 


2 


1 


1 


1 


5 


5 


5 


4 . 


8 


6 


3 


5 


4o 


300 


347 


% 




1 


1 


4 


1 


4 


1 


4 




2 


1 


5 


3 


5 


7 


4 


T 


5 


2 


12 


130 


1610 


172O: 


4 


3 


< 1 


2 


5 


9 


2 


4 


Krj 




1 


3 


6 


9 


10 


51 


4 


24 


17 


5 


24 


4-50 


84oo 


12400 


4 




4 


< 1 


1 


1 


l 


4 


4 




5 


1 


2 


5 


27 


42 



(Table I continued on next page ) 
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FUNCTION 


X = 


X = -3 


x = .6 


X = >9 


X = .95 


x = ,98 


x s= .985 


*8 


10 


8 


25 


50 


1500 


37000 


16200 


4 


< 1 


2 


1 


it 


1 


3 


2 


K 9 


16 


8 


7 


5 


2k 


133 


c620 


*9 


87 


57 


50 


105 


if200 


42000 


300,000 


^0 




1 


1 


2 


2 


1 


4 


K 10 




50 


50 


25Q 


7& 


*&-■■ 


11QQ 


^0 


2 


120 


177 


270 


10,100 


1,1^0,000 


■worthless 



















LOCATION 



7 
8 

9 
.10 

11 

12 

13 
14 

15 
16 



ORDER 



40 3F 

41 5F 
iH F 
X5 F 

42 6jL 
10 20F 
42 SF 

49 S3 
L5 3F 
10 IF 
40 1S3 
F5 F 
40 4F 
L5 3L 
42 13L 
42 28L 
42 62L 
L5 13L 
42 1?L 
L4 4f 
42 13L 
F4 4f ^ 
42 19L 
F5 4f 
40 4f 
FO 5F 
32 20L 

50 ( )F 
7 J 3F 
40 IF 
L5 4f 
FO F 

50 F 
66 4f 



NOTES 



PAGE 1 



V 6 



Form 20 



By 7 and 
10 



Store argument in 3 

clear k storage location 

n ° 

Create a zero 
Plant link 



_ Pl^nt. k 



n 



Set 1/2 P Q = 1/2 in S3 

- Set 1/2 P x = 1/2 x in 1S3 

1 Set k = 1 

Set initial 
addresses to S3 



Set address of P. 



k-2 



Set address of P. 



k-1 



Set address of P, 



Advance k -» k .+ 1 



k - k 
n 

@ k > k , Done with P r s 
g) Not done with P's 
x 1/2 P^ -> 1 



k-1 



in Q 



LOCATION 


ORDER 


NOTES 


PAGE 2 


IT 


L5 IF 










LO ( )F 


By 9 


complete l/2 P, 




18 


40 2F 
7J 2F 








19 


1,4 IF 










40 ( )F 


By 11 


Store l/2 P. 
1 k 




20 


22 8L 




Repeat cycle 






kl kY 


From 13 


Set k counter to 




21 


L5 28L 
Lk kY 


From 32 


Advance 




22 


k2 28L 
k-2 29L 




addresses 




23 


iH F 
50 F 




Clear 
Clear Q 




IX 


L5 ^F 
Yk kY 




J-V- 2 - inA 




25 


00 3^F 
66 70L 








26 


S5 F 

50 26l 










r(2k+D , n P 


27 


26 71L 
19 IF 




-1 -2 

2 in A 




28 


50 2F 
74 ( )F 


By 7 and 








\ 


2k + 1 1ti q 


Htr 


29 


00 IF 


22 


\ — — »» 






kO ( )F 


By 22 


Store K 




30 


F5 kY 
kO kY 




— advance k counter 




31 


FO 5F 
32 32L 




k - k 
n 

Done normalizing 




32 


26 21L 
50 3F 


From 31 

i 


@ Continues 





V 6 



LOCATIOH 


ORDER 


NOTES PAGE 3 


33 


71 3F 




2 . A 
-x in A 




L4 691 




2 . 
1 - x in A 


34 


40 IF 
S5 F 




NSP (1 - x 2 ) -> 1 


35 


40 F 




LSP (1 - x 2 ) -» 




50 35L 






36 


/TTf J J. J-l 

40 8f 






Jl - x in 8 


37 


4l F 




Clear 


38 


F5 F 
4o 4f 




-Set k = 1 




40 7F 


From 68 


Set m = 1 


39 


F5 62L 








42 42L 


From 65 


Plant address of KT 


4o 


Fk F 




Plant address of K~ 




42 62L 






41 


£0 4f 
42 57L 




Plant address of K=V, 


42 


50 3F 








7J ( ' )F 


By 39 


-xK?" to 6 


^3 


4o 6f 








L5 4f 






44 


L4 7F 
4o 9F 




- (k + m) -*9 


^5 


FO F 
40 IP 




-(k + * -l)-*l 


46 


l4 4f 








LO JF 




-( 2k - 1) -»Q 


47 


10 39? 




_i 




75 9F 




"_ (k + m )(2k - 1) ->2 


48 


S5 F 






r 


40 2F 


_, 



V 6 



LOCATION 


ORDER 


NOTES PAGE k 


^9 


F5 ^F 








iM- 4F 




- (2k + 1) -»Q 


50 


10 59F 








75 IF 




"_ lA (k + m - 1) (2k + 1) inA 


51 


00 37F 






52 


66 2F 
S5 F 
50 52L 








- 1/2 


(k+m-l)(2k+l) , n 
(1+m) (2k -1) -* 10 


53 


26 71L 








kO 10F 






5^ 


F5 4f 








LO 7F 




(k - m + 1) in A 


55 


50 F 




Clear Q 




66 9F 






56 


S5 F 

50 56L 




~P 




- m + 1 . 

= in 2 

k + m 


57 


26 TIL 




wl 


r." 


50 ( )F 


By hi 


_ Ci *« 


58 


75 10F 
00 IF 






r(k+m-l)(2k+l) . A 
— I?* — — — *•*£ in A 


59 


kO 10F 




J y (k+m)(2k-l) 

HSP [ " '] in 10 




S5 F 




LSP [ " ] in A 


6o " 


50 2F 

7$ 6f 






6l 


lA 10F 

66 8f 






62 


S5 F 






63 


ko'( )F 
^5 7F 


By 8 and 


Store K? 




kO 7F 




— Advance m -» m + 1 


6k 


FO 4F 




m - k in A 




56 66l 




© m > k, done this k 



V 6 



LOCATION 


ORDER 


NOTES PAGE 5 


65 


•L5 62L 




O Not done this k 




22 3% 




Repeat 


66 


F5 kF 


From 6^- 






kO kF 




- advance k -» k + 1 


67 


FO 5F 




k - k in A 
n 




32 ( )F 


By 2 


Q Done, enter Link 


68 


F5 F 




® not done j 




22 38L 


_ 


Set m. - 1 and. repeat 


69 


80 F 


_ 






00 F 


- 


-1 = 1 


70 


kO F 








00 285 398 163 397J 


8x ■& 
-t 32 " T~ 


71 


50 72L 








26 999F 




Interlude 


72 


00 F 




to load 




00 71L 




m ^ 


73 


00 F 




Routine R - 1 116 




26 71L 




at 7^L 


26 IN 






Routine El - 116 






Square Root 

1 







V 6 



